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We present a two-loop computation of the beta functions and the anoma- 
lous dimensions of a 75- Yukawa model using differential renormalization. The 
calculation is carried out in coordinate space without modifying the space- 
time dimension and no ad-hoc prescription for 75 is needed. It is shown that 
this procedure is specially suited for theories involving 75, and it should be 
considered in analyzing chiral gauge theories. 
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I. Introduction 

The differential renormalization (DR) method |1| has proven to be a sim- 
ple and successful way to treat ultraviolet divergences in a renormalizable 
quantum field theory @ , || , , [(J • The method is developed in coordinate 
space and does not change the dimensionality of space-time. It consists of 
two steps. In the first one bare amplitudes, which are too singular at short 
distances to admit a Fourier transform, are replaced with derivatives of less 
singular functions. In the second one, derivatives are integrated by parts. In 
doing the latter, ill-defined surface integrals arise: they correspond to coun- 
terterms that guarantee finite (renormalized) Green functions and unitarity 
[|]. Thus, by discarding surface integrals, one simultaneously regularizes and 
renormalizes divergent amplitudes. 

One of the motivations to develop DR was to find a consistent way of 
treating theories with dimension-specific objects [[J. Theories involving 75 
are of this type and are going to be the subject of this letter. 

Dimensional regularization is the most practical and better understood 
regularization method used in quantum field theory but it runs into diffi- 
culties with theories involving 75. There is a large amount of literature and 
controversy (see p] for a review) over the correct prescription for 75 in dimen- 
sional regularization. There are basically three different ways to deal with 
75 when dimensional regularization of the Feynman integrals is employed, 
each one of them giving rise to a different regularization scheme. The first 
one is called naive dimensional regularization || and uses an anticommut- 
ing 75 and the 4-dimensional Dirac algebra. The second one goes under the 
name of dimensional reduction M, where both the Dirac algebra and the 
vector fields are regarded 4-dimensional and split into D-dimensional bits 
plus 4-D scalars. The third one is the 't Hooft-Veltman proposal |10| and 



uses a non-anticommuting 75 in D dimensions. It has been shown that the 
only systematic, uniquely fixed and consistent regularization procedure is 



the original prescription of 't Hooft and Veltman |ll|,[12|,|7fl- The other 



two schemes present manifestly algebraic inconsistencies; they may work for 
some low-loop calculations, specially when there are no closed odd parity 
fermion loops in the theory, but the validity of the computations to all loop 
orders is not guaranteed. However, in the 't Hooft-Veltman prescription 
the computational simplicity of dimensional regularization is lost and (D- 
4)-dimensional or evanescent counterterms have to be taken into account to 
obtain the correct results. This scheme is rarely used due to these computa- 
tional difficulties. 

It is our aim to show that DR provides a convenient way to handle theories 
with 75. Since in DR the dimensionality of space-time is kept unchanged, 
the 4-dimensional Dirac algebra and the standard 75 definition are main- 
tained. As a result, the algebraic consistency is ensured. The simplicity of 
the computations is also another salient feature of DR. 

As an example, we will study a massless Yukawa model with lagrangian 

C = -(dip) 2 - —(p A + ii/j-fd^ + ie(filj^. (1) 

This model has been analysed in the literature using dimensional regular- 
ization with two of the above prescriptions, namely, the naive dimensional 
regularization ]T3 , and the 't Hooft-Veltman prescription [I3|. Here we use 



DR to compute the beta functions and the anomalous dimensions up to two 
loops thus avoiding the mentioned difficulties inherent to dimensional regu- 
larization. In sections II and III some details of the calculations are carefully 
explained. Section IV shows the results and section V is devoted to the 
discussion. 

II. One-loop order 

We will start by studying the model (|]) at one-loop. This will give us the 
chance to recall the basic techniques used in DR. 

We will work in Euclidean space, the Euclidean Feynman rules and some 



other conventions being given in Figure 1. Figure 2 shows the 1PI divergent 
graphs of the model at one and two loops. 

The bare one- loop contribution to the fermion self-energy (Fig. 2.1) is 
given by 

lj P \ 2 1 1 a 1 

*!-<*. y) = - Uy^T^W 1 - Wlf = -&*&=&• (2) 

where (let us emphasize) the standard 4-dimensional Dirac algebra has been 
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used, and g = jj^y. From now on, we will use the convention that derivative 
operators act on the first variable of the fraction, so that di(x — y)~ 2 = 
—di(y — x)~ 2 . Clearly, the singularity of the bare amplitude (^) lies in the 
l/(x — y) 4 factor, which diverges too strongly when x — > y so as to have a 
Fourier transform. Let us consider for a moment x ^ y. Then, the equation 

I 

1 1 lnQr -y) 2 M 2 

{x-yY 4 (x-y) 2 ' U 

is an identity. Now comes the key point: DR defines the renormalized ampli- 
tude as the result of using Eq. (|]) (also at x — y) with the proviso that, when 
computing the amplitude in momentum space, the derivatives acting on the 
r.h.s. above are formally integrated by parts and the additional prescription 
that surface integrals are discarded. The constant M has dimensions of mass 
and plays the role of a subtraction point. More explicitly, Eq. (|3p gives for 
the renormalized one-loop contribution to the fermion self-energy: 

g Hx-yfM 2 
8n z (x — y) z 

Now we can compute the amplitude in momentum space (setting y=0 due 
to translational invariance) 

i«./^(>^). (5, 



Following DR, we integrate by parts discarding the surface term and are left 
with 

«») = ->^' (6) 

where 7 = 1.781072... is the Euler constant. 

The bare one-loop boson self-energy (Fig. 2.2) is 

A!"(x, V) - "U* {^^-2-^^} = ^^. (7) 

The singular factor l/(x — y) 6 in (0) now produces, upon Fourier transform, 
a quadratic divergence. To cure it, we have to extract two derivatives more 
than in @. To this end, we use (see Appendix A of ref. [|TJ) 

1 1 \n(x-y) 2 M 2 

= □□— - — , (8) 

(x -yf 32 {x-y) 2 ' l } 

thus obtaining the renormalized amplitude 

a / 9 \n(x-y) 2 M 2 

We could compute now the renormalized amplitude in momentum space, as 
in the previous case. 

We next renormalize the one-loop Yukawa vertex (Fig. 2.3). The bare 
amplitude is 

V*™{x } y,z) = 7^75^7-^75^7-^75- ' 



(47r 2 ) 3 ' (y — x) 2 " {x — z) 2 (y — z) 
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= -ie^7bF(x,y,z). (10) 

To regulate this graph we notice that the function F has a singular region 
when x ~ y ~ z. Outside this region we can integrate by parts and obtain 

n, d C ! 9 11 

■F(s, 2/, ^) = 7i7j 
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<9x* [ (x — y) 2 (y — z) 2 dxi (x — z) 

.,.., I * L_ (11) 



The first term on the r.h.s. has a well-defined Fourier transform, as can 
be seen using dimensional analysis. In turn, the second term is singular at 
x ~ y ~ z\ notice that it contains the laplacian Y^didj and that 

n-^-- = -4n^(x-z), (12) 

(x — z) 2 

so the second term in ([11]) can be renormalized using ([!]). Finally we can 
write the DR renormalized graph as 

V( \ *l i d_( 1 d I \ 

v l{ x,y,z) - 47r4 75|7i7i 9a;i [( x - y )2( y _ z y dx j ( x _ z y ) 

- nWHx-z ^-^f ). (13) 

[x-yY J 

The 4-point boson function has two contributions at one-loop: the bubble 
(Fig. 2.4) and the box (Fig. 2.5). The bare bubble is 



t*?(x,y,z,w) ' 



32tt 4 



8^'ix — y)5^ 4 '(z — w)- — + permutations 

(x — Zf 

(14) 
and is easily DR renormalized using (H): 



Ti a (x,y,z,w) = — 2h 2 5 ( - 4: \x — y)8^\z — w)0 — — \- permutations, 

[x — zy 

(15) 



where h 



A 



16tt 2- 

The bare box 
rfr°(x,y,z,w) 



(ieY 



(A,r tr {^j^^^j^^^j^^^V^wA +perm - 

(16) 
is a bit more complicated due to the presence of indices but the DR program 
can be straightforwardly carried out using the Dirac algebra and identifying 
the singular factors. In this case, we integrate by parts twice and use (||) to 
obtain the following DR renormalized amplitude 
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T lb (x,y,z,w) = 

tr{j a j b j c j d }—- -d b - -d c - —d d 



7r 4 """'"' ' Qya y^y _ x y (y _ w y ^ w _ z y (z — x) 2 ) 

it 2 {y — x) 2 dz d \(z — x) 2 dz d (z — w) 2 J 

+ I6g 2 5^\y - w)6^(z - w)D ln( f - W \f + perm. (17) 

(x — w) z 

This analysis shows that all one-loop UV divergences are cured by first 
isolating the singularities and then by using Eqs. (|3]) and (|8|), and the inte- 
gration by parts prescription. In this process an arbitrary parameter with 
dimensions of mass enters in a natural way. In writing the different renor- 
malized amplitudes, one has no a priori reason to use the same M for all 
graphs but we will use a subtraction scheme where all the M's are equal. 
This can be done since there are no Ward identities in the theory that force 
us to introduce different mass scales [I5|. This parameter has the physical 



meaning of a renormalization scale, as can be explicitly seen when the renor- 
malization group equations are tested on the renormalized amplitudes (see 
Section IV). 

III. Two-loop order 

For the sake of brevity, we will not discuss here all the graphs of the two- 
loop order but we will present an overview of the computations, which involve 
almost the same techniques used in the one-loop order, that is, isolating the 
singularities of bare amplitudes and using DR identities. 

Two-loop bare amplitudes present two singular regions. The regulariza- 
tion of these amplitudes is performed from the sub divergence, treating it in 
the same way that it was done at one-loop, to the overall divergence. Overall 



divergences require new DR identities. In our case, we only need 

lnx 2 M 2 1 \n 2 x 2 M 2 + 2\nx 2 M 2 

— □- 



x 4 8 x 2 

lnx 2 M 2 1 ln 2 x 2 M 2 + 51nx 2 M 2 , . 

= an . (18) 

x e 64 x 2 V } 

In this way two-loop 1PI are renormalized. As an example, we examine one 
of the diagrams which gives contribution to the fermion self-energy (Fig. 2.6), 
with bare amplitude 

S h 2 Z e {x,y) = 

e 4 1 f d A ud A v \ 



(47r 2 ) 5 " (x — y) 2 J (y — v) 2 (x — u) 2 \ (v — u) 
Ag 2 ^ 1 r d A ud 4 v 




7r 6 (x — y) 2 J (y — v) 2 (x — u) 2 (u — v) & 

The first singular region of this amplitude is the corresponding to the sub- 
divergence, when u ~ v, and is cured in the same way that it was done at 
the one-loop order (|9|). The laplacians in front of ln(ti — v) 2 /(u — v) 2 are 
then integrated by parts, discarding the surface terms. Then we can perform 
easily the integrals over the internal points. The second singular region cor- 
responds to the overall divergence, when x ~ y, and is cured using (|T^) to 
finally obtain the DR renormalized amplitude 

Q / s 9 2 a ^ ln 2 (x-y) 2 M 2 + 31n(x- j /) 2 M 2 ^ 

s*(*. v) = -^ [ D ^z^ J • ( 2 o) 

Amplitudes with well located sub divergences are cured in a straightfor- 
ward way, so we will not discuss them. There are other diagrams which re- 
quire a supplementary effort: those that contain an overlapping divergence. 
In those diagrams it is impossible to separate out the two singular regions. 
Nevertheless, overlapping divergences do not represent a serious difficulty in 
real-space computations, because the external points of the amplitude can be 
kept separated until the regularization of the sub divergences is accomplished. 



Let us take as an example of overlapping divergence the amplitude cor- 
responding to Fig. 2.10, 

A*T(x,y) = -j^str{ lalblcld }x 

du A dv 4 \ d a - -d b - -d c - -d d - - \ . (21) 

[ \x — u) [u — y) z [y — v ) z (v — x) z (v — u) z J 

In this amplitude the singular regions are u ~ v ~ x and u ~ v ~ y. To 
renormalize it we first extract total derivatives over the external points. We 
are then left with three terms. We notice that in two of these terms there 
is a laplacian acting on a fraction that allows us to obtain a delta function 
and perform one of the integrals over an internal point. These two terms 
present singular regions when u ~ x in one case, and when u ~ y in the 
other, and they are renormalized using (H) first, and (||) and (pl|) afterwards. 
Also the third term can easily be renormalized. After some algebra, the DR 
renormalized amplitude is expressed as 



g 2 ln 2 (x-y) 2 M 2 + \n(x-y) 2 M 2 
27r 2 (x - y)' 



A 2b (x,y) = ^ ' y> si T~ y ' • (22) 



The very same procedure is used to renormalize overlapping divergences 
of the Yukawa vertex and of the 4-point boson function : extract total deriva- 
tives and use the properties of Dirac matrices during the computation. In 
these diagrams, and due to the presence of indices, there are some integrals 
that are difficult to solve analytically, but fortunately they are finite, and they 
do not need to be renormalized. For our purposes, to find beta functions and 
anomalous dimensions, we do not need them in a explicit form. 

As an example of overlapping divergence of a vertex, we next examine 
the most difficult two-loop Yukawa vertex (Fig. 2.16), with bare amplitude 

V 2 b r(x, y, z) = ___ 757i7i7ryifc x (23) 



d A ud A v I di — dj — d k — d t — — 

[ [x — u) z [u — y) z [z — v) A [x — v) z [u — z) z (y — vy 



To proceed, we first integrate by parts over d/dx l . Power counting shows 
that only the integrated part is too singular so as to have Fourier transform, 
so we foan our attention only on that part. The laplacian acting on l/(x — u) 2 
gives a delta function, and we can then perform the integral over u. The final 
integral to be analysed is 

47T 2 7 i 7 fe - T-rdj- -r / (fv- t~t- r=8 k 



x — z) 2 (y — x) 2 J (x — v) 2 (y — v) 2 (z — v ] 
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Now d/dyi is integrated by parts. Once again, only the integrated part 
needs renormalization. This last one is separated into traceless and trace 
parts. The traceless combination of derivatives is finite, and the trace part is 
straightforwardly DR renormalized. After some algebra, the final result can 
be expressed as 
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V 2e {x,y,z) = y^75S -7j7i7/7fc x 



d /„ 1 1 1 1 



d 4 ud 4 v^—r ( di- — ^-^dj- — ^-^d, 
) 2 J {u-y) 2 

d ( 1 d 



dx l \ l (x — u) 2 J {u — y) 2 (z — v ) 2 (u — z) 2 {y — v ) 2 (x — v ) 2 
+ 4vr2 7j7fc7T- ( 7- -^7-, — ~^1TT K ( Z -x,y-x) 



' dyi \{x — z) 2 (y — x) 2 dz k 
4 2 1 / d 2 Sjt d 2 \ _ ,. 

3 (x — z) 2 (y — x) 2 \dyidz k 4 dy-dzj 



+ 8tt 4 



1 ln(x - y) 2 M 2 ln(x - z) 2 M 2 
16 (x — y) 2 (x — z) 



2 



+ 1 d l \n{x-y) 2 M 2 d 1 

4 <9x^ I (x — y) 2 cfe^ (x — y) 2 (x — z) 2 

tt 2 (4) Am ln 2 (x-y) 2 M 2 + 21n(x-y) 2 M 2 
+ y^x-^D j— — + (?/<-*) 

Although the final expression seems to be somewhat complicated, specially 
for taking its Fourier transform, the computation of Md/dM on it, which 
is what we need to find the beta functions and the anomalous dimensions, 
presents no difficulty. 

The remaining overlapping divergences of the Yukawa vertex and of the 
4-point boson function can be treated in a similar way, and we will not 
explicitly discuss them here. 



IV. Renormalization group constants 

After the computation of the renormalized amplitudes, it is easy to find the 
values of Md/dM acting on all the 1PI graphs. We present them in Fig. 
3 in a pictorial form. With those values, it is possible to check that the 
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renormalized amplitudes satisfy renormalization group equations 

[ M Mf + ^T g + ph M - n ^ ~ n ^) r(w) = °' (26) 

When we substitute the values of Md/dM in fl26|), we find the following 
values of beta functions and anomalous dimensions 

a, = io 9 2 + i^-4V-f«? 3 

A = 3,r + 8,, 9 - 48^-^-12/^ + 28ft/ + 38 V 

1 1 13 

lv = 2g + —h 2 -5g 2 ; ^ = -g-—g 2 . (27) 

These results coincide with the ones given in the literature, where the same 



computations were carried out using naive dimensional regularization [[13 
and the 't Hooft-Veltman scheme Hl^lFJ 

This model does not present closed odd parity fermion loops, so one does 
not expect the naive dimensional regularization to fail here. Nevertheless, 
one would prefer to use a procedure which does not break the algebraic 
consistency of the computations and such that one could always rely on its 



validity. These problems do not arise in [ 1 1 1 , but there the cumbersomeness 
of the computations, where evanescent counterterms had to be taken into 
account, shows that the 't Hooft- Veltman scheme is almost impractical. 

We must remark, once again, that computations using DR are simple and 
algebraic consistent. 

V. Discussion 

We have applied DR to compute the beta functions and anomalous dimen- 
sions up to two loops of the Yukawa model defined in Eq. ([!]). 

1 There is a 2 factor between the beta functions obtained in |14j and our results, due 
to the use of a different convention. 
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DR presents some salient features. DR is just a prescription to extend 
distributions which are ill-defined in some points to distributions well-defined 
in all space-time. Bare and DR renormalized amplitudes coincide for non- 
singular points, so this procedure represents a minimal change in the theory. 

DR is a real-space method, and computations of multiloop integrals can 
be done with ease. Let us recall here that the "method of uniqueness" 



n^)[r71)lI8H which also works with coordinate amplitudes but uses dimen- 
sional regularization, has already demonstrated the powerful calculational 
possibilities of space-time regularization methods. DR also shares that com- 
putational simplicity. The method is specially suited to obtain renormaliza- 
tion group functions and renormalized amplitudes in configuration space at 
high loop order. The computation of these amplitudes in momentum space 
requires to Fourier transform and perform some difficult but well-defined 
integrals. 

In theories in which 75 is present, DR preserves the algebraic consistency 
of the computations, since it does not alter the dimensionality of space-time, 
and the standard 4-dimensional Dirac algebra can be used. 

It is worthwhile to mention here that DR has been used in gauge theories 
[||], [|J, and that it reproduces correctly the chiral anomaly ||. In gauge 
theories DR enforces the introduction of several mass scales, and in order to 
preserve the Ward identities certain correlations between them are required. 
The anomaly results when the Ward identities overconstraint the value of 
the mass parameters. 

There is still no proof of the validity of DR to all orders of perturbation 



theory. However, it has very recently been shown [15] that there is a relation- 
ship between dimensional and differential regularization and renormalization 
in low loop graphs, and one expects that one could find a systematic appli- 
cation of these relations at higher order graphs and then yield a consistent 
proof of DR. One difference between the two procedures is that DR does 
not modify the space-time dimension, and then it allows for a natural and 
consistent treatment of dimension-specific theories. 

13 



DR seems to be a suited and natural method to treat theories with 75 
and, with no doubt, deserves further consideration to evaluate chiral gauge 
theories. 
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Figure Captions 

Figure 1. Euclidean Feynman rules and other conventions used in the com- 
putations. 

Figure 2. One and two loop 1PI diagrams. In parenthesis, the number of 
diagrams of that type. 
Figure 3. The value of Md/dM acting in all the considered 1PI. 
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